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Abstract 

Nontopological solitons can be formed during a phase transition in the 
early universe as long as some net charge can be trapped in regions of 
false vacuum. It has been previously suggested that a particle-antiparticle 
asymmetry would provide a source for such trapped charge. We point out 
that, for the model and parameters considered, statistical fluctuations 
provide a much larger concentration of charge, and are therefore, the 
dominant source of charge fluctuations in solitogenesis. 
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Nontopological soliton solutions in classical field theories have appeared in many 
forms since they were first introduced by Rosen 1 and by Friedberg, Lee, and Sirlin. 2 
Examples include Q balls, 3 quark nuggets, 4 cosmic neutrino balls, 8 and soliton stars.® 
The simplest nontopological soliton (NTS) solution involves a real scalar field, tr, 
and a complex scalar field, <f>, with Lagrangian 7 

U{\<t>\,<r) = — cr*) 3 -h — <r 0 ) 3 H- y (<r - co) 3 ^ + A, (1) 

where the constant A is adjusted to give 17(0, tr) = 0 at the global minimum of the 
potential. The classical potential for a has two local minima. At the global miminum 
(<r = <r_ = —[(1 + 2Aj/Ai)/2 + [(1 + 2A 2 /Ai) 2 + 8A 2 /Ai] 1 / a /2]<r 0 ) the <f> field has a 
mass = h(<r_ — <r 0 ) 2 , while at the local minimum (<r = <r 0 ) the field <p is massless. 
The nontopological solition solution describes a configuration of massless <f> particles 
trapped inside a region with a = <To, separated from the true vacuum a = <r_ by a 
wall of thickness ~ cr^ 1 . The energy of the NTS configuration is® 

+ + O(R 2 \\ /2 <r 0 ), (2) 

XL O 

where Q = \N+ — N^\ is the “charge” contained in the spherical NTS of radius R, and 
A is given by A = -Ai(<rl - o-q) 2 / 8 - A 2 (<r_ - <ro)V 0 /3. The three terms in Eq.(2) 
represent the kinetic energy of the confined massless <f> field, the false vacuum energy 
of the NTS, and the surface energy of the wall separating the interior NTS region 
from the true vacuum. We will ignore the wall contribution in our analysis. 

Minimizing the energy of the NTS configuration results in an NTS of mass and 
radius M = (4:Tr/3)y/2Q 3 t A A 1 t* and R = (Q/ 4A) 1 / 4 . This mass should be compared 
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to the mass of Q free (f>'s in the true vacuum, Mfr ee = Qm $ = Qh}t 2 \<r_ — a 0 \. The 
NTS configuration will have a lower mass, and hence be stable, for charge Q greater 
than some minimum charge, given by 

_ 1231 A , 0 . 

Qum = ~h 5 >- -*.)*' (3) 

In this paper we will study in detail the case Aj/Ai = 0.15. For this choice of A2/A1, 
A = O.6Ai00 , Qmin = 18Ai/h 2 , and A/min = 46(Ai/A 3 / J )cro. 

A scenario for the cosmological origin of NTS was proposed by Freeman, Gelmini, 
Gleiser, and Kolb 8 (hereafter, FGGK). In the FGGK scenario, there is a critical 
temperature, Tc — 2<r 0 , below which the Universe divides into domains of true (<r = 
cr_) and false (<r = <ro) vacuum. The characteristic size of these domains is determined 
by the correlation length, £, of the cr field at the transition. At high temperatures 
thermal fluctuations can cause a correlation volume to make the transition between 
the two minima. These fluctuations freeze out at the “Ginzburg” temperature, To- 
FGGK estimate To by the criterion that To is equal the maximum free energy of the 
correlation volume in the transition Fm = UmV £ ( Um is the maximum value of the 
potential in the region < a < (To). For A2/A1 = 0.15, To = 1.3<ro/Aj^ 2 . Of course 
To can never be larger than Tc — 2<r<). 

At To , the probabilities of being in the false vacuum, p(<7 0 ), and true vacuum, 
p(<r_), are Boltzmann distributed according to the difference in free energies of a 
correlation volume in the different minima 

^ = expt-Af/iy = exp[-AVr 0 l (4) 

(recall that U{(rS) = 0 by the addition of A). If p(o r o)/p( £r -) ^ 0*3 then only finite 
regions of “false” vacuum will be populated. If the regions of false vacuum contain a 
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net charge Q = \N <t> - Nj | > Qmin, true false vacuum regions can be stabilized and 
evolve to nontopological solitons. 

The probability that a false vacuum region contains a charge Q > Qmjn is the 
subject of this paper. 


II. CHARGE FLUCTUATIONS 


FGGK assumed that the net charge in a region was proportional to a cosmic asymme- 
try, like baryon number, between <}> and <j>. This cosmic asymmetry can be expressed 
in terms of an asymmetry parameter 77: 


I n * ~ n *\ 

V = — — — • 

n* + n ? 


(5) 


In this paper we demonstrate that if 77 < 0.5, Poisson fluctuations will dominate the 
probability distribution, and the number density of NTS’s produced will be indepen- 
dent of 77, even in the limit 77 — ► 0. 

We will denote by N4, the mean number of <p's in some volume. The proba- 
bility of finding the actual number N# of </»’s is Poisson distributed: N#) = 

In the limit of large Nj,, the distribution will be Gaussian, with mean 
and variance /x = a 3 = N+: 


P(Nf, N t ) = — = «P [-(N t - N t y/ 2 Ni\ . (6) 

We will later discuss the validity of the Gaussian approximation. An expression 
similar to Eq.(6) obtains for the probability of finding a number Nj of 4 >' s if the mean 
is Nj. The total number, N, and charge, Q, defined as 


N = N 4 + JVj Q=\N+- N ; | 


(7) 
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will also be Gaussian distributed; with means N = N# + Nj, Q = | N# — Nj\ = t)N] 
and variance cr 2 = N. Therefore, the probability of finding a charge Q in a volume 
containing a mean number N of (<f> + 0)’s is 9 

P(Q, V) = ^=^exp [-((?- ,57)7277]. (8) 

As described by FGGK, below Tg the Universe divides into cells of correlation 
volume V( ~ (47r/3)£ 3 . Adjacent cells of false vacuum form “clusters” with density 
per unit cluster of 

f(r) = 6r -1,5 e -cr (9) 


for volume V = rV, j. The constants b and c are unknown. Scaling arguments imply 
that c — ► 0 as p(<Tq) — * p e (where p e is the critical probability for percolation, p e ~ 1/3) 
and b — ► 0 as p(<r 0 ) — > 0. It is expected that b and c are of order unity otherwise. The 
number density of r-clusters produced in the transition is simply n(r) = /(r)V^ _1 . In 
a volume V = rU(, the mean number of ( <f> + <f>y s is N = rN 4, where N( is the mean 
number of ( <f> + <f >)' s in a correlation volume 

The number density of false-vacuum domains with charge Q is simply given by 
nq = S^l 1 n(r)P(Q; N = rN {), where n(r) = /(r)V^ 1 as before, with /(r) given by 
Eq.(9). Approximating the sum over r by an integral, 10 nq becomes 

V^nq = 7 dr r ~ 2 exp \-( v 2 N ( /2 + c)r - Q'fiNfr] 

= + c) 1/3 *i [( V2Q/N\ ,2 )( V 2 N</2 + c) 1 ' 2 ] , (10) 

where Ki(z) is a modified Bessel function of the second kind of order one. For large 
argument, the expansion ■K’i(z) — ► e~ z \Jir/2z gives 10 
TV 1 / 3 

V i n Q = 6 Q3/2^ 3 + 2c /^) 1/4 ex P [Qv - QW + 2c/77{) 1/2 ] . (11) 
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The expression for V^nq can easily be converted into the ratio of the number 
density of <?’s to the entropy density, s = 2tt 2 g,T 3 jAb. Using £ = A ^Tq 1 , the 
correlation volume is V( = 47t£ 3 /3 = 47t/3A 3 Tq. Assuming the <f >' s are relativistic at 
T = Tq, n-^ = C(3)Tg/7t 2 , and at T = Tq 


N t = ( n * + ~ 3 - 


( 12 ) 


Since = 87r 3 <7./135A 3 a, 

tiq 135A 3 6 


y 0 = ^ = 


8((3) 


!■/> 


37rA 3 c\ 1/,< 


(i , 3 ? rA i c \ 

r 4C(3)j 


3 87 r 3 g+Q 3 / 2 [ 37rAf J 

QV~Q (v 2 + 37tAjc/4£(3)) 1/2 


x exp 


= ^^(, I + 1.96M)'/*exp[ ( ?,-<?(,> + 1.96Ajc)V i ]. (13) 

There are two interesting limits of Eq.(13). In the limits tj 2 1.96A 3 c, and 

tj 2 1.96A 3 c, Yq becomes 


Y _ { (0MbXl /4 /g,Q^ 2 )exp(-lAQXl /:t c^ 2 ) t) 2 < 1.96A 3 c , , 

Q \ (0.54b\^ 2 r] 1 / 2 /g t Q 3 / 2 )exp(—0.98\lQcT}~ 1 ) tj 2 » 1.96A 3 c ' 

Since Yq decreases exponentially with Q, the most abundant NTS will be the one 
with the smallest allowed charge, Q = Qmin- 

Note that in the “large” tj limit, we essentially recover the results of FGGK. How- 
ever, this case is only relevant for tj larger than of order unity. A much more likely pos- 
sibility is that the “small” tj limit is the relevant one, and that Yq ~ 10 - 3 QMiN a e -9MIN . 


III. NUMERICAL RESULTS AND CONCLUSIONS 

In the previous section, three approximations were used: 1) Gaussian rather than 
Poisson statistics, 2) the sum over r-clusters was replaced by an integral, and 3) 
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the large-z expansion of the Bessel function K\(z) was used. In this conclusion 
section we present some numerical result and discuss the range of validity of the 
above approximations. 

Clearly for “large” Q, Q > 10-20, Gaussian statistics will be a good approxima- 
tion. In Fig. 1 we compare an integration over r of Gaussian statistics, Eq.(10), to the 
more accurate sum over r of Poisson statistics. The Gaussian results are presented 
for rj =0, 0.25 and 0.5, while the Poisson results are given for rj = 0 only. It is clear 
that the Gaussian approximation is an adequate one. Integration over r rather than 
summing also introduces only a small error. 

In Fig. 2 we present the large-z expansion of the Bessel function in Eq.(ll). Com- 
parison of Fig. 1 and Fig. 2 shows that for ^/2QNJ 1 ^(tj 2 N^/2 + c) 1 / 2 > 2, the 
expansion is accurate. In Fig. 2 we also show for comparison the results of FGGK for 
n Q V v Clearly it is a serious underestimate for tiq unless r j a 1.96A 3 c. 

We conclude by illustrating the importance of the calculation of Yq. We use the 
example discussed in the introduction, Aj = 0.15Ai, which gives Qmin = 18Ai/h 2 , and 
•^(Qmin) = 46Ai<r 0 //i 3 / 2 = 2.5<2MiNh 1 / 2 <ro* Assuming that the contribution to 17 from 
NTSs is dominated by those with Q = Qmini Ynts — 10 -3 Qmin e - ^ MW , the present 
NTS energy density is pnts = YntsMntsSq, where is the present entropy density, 
so =2800 cm -3 . Comparison of pnts to the critical density, /)c = l*88xl0 _29 h 2 g cm -3 , 
where ho is the Hubble constant in units of 100 km s -1 Mpc -1 , gives 

^NTshl = 10 6 /t 1/2 ( -fTTr ) ~Tj5~ • ( 15 ) 

Vmin 

For NTSs to be dynamically relevant today, fl^rsh 2 should be in the range 10" 2 < 
ft NTS^l < 1. Relevant values of or equivalently Ai/h 2 , are shown in Table I. 

The conclusion of this paper is that statistical fluctuations are the dominant source 
of charge fluctuations in solitogenesis, not a cosmic asymmetry as assumed by FGGK. 


6 



* 


The resulting Qnts is independent of 77, so long as rj < 0.1. Finally, reasonable values 
of Qmin give SInts in «• dynamically interesting range. 
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TABLE I 



0-0 = 10 19 GeV 

<tq =1 GeV 

ftNTsh 0 

<?MIN 

Ai/h 2 

Qmin 

Ai /h’ 

10" 2 

60 

3.3 

17 

0.94 

1 

55.5 

3.1 

12.5 

0.7 


r 
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FIGURE CAPTIONS 


Figure 1: A comparison of Poisson (Footnote 9) and Gaussian probabilities (Eq. 10) 
as a function of Q. N^ = b — c=l was assumed. 

Figure 2: The result of the large -2 expansion of the Bessel function in Eq. 11 is shown 
by the points marked This Work. Comparison of these points with the corresponding 
points in Fig. 1 shows that the large -2 expansion is a good approximation. Also indi- 
cated by the points marked FGGK are the results of FGGK 8 which ignored statistical 
fluctuations. 
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